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摘要 : 光学 窗 波 段 观 测 能 够 表征 发 射线 天 体 清 晰 的 结构 ， 是 研究 发 射 
x. BAHU 区、 星际 介质 等 方面 的 研究 具 
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度 来 调整 图 像 对 齐 ， 然 后 进行 两 幅 图 像 的 相 减 ， 但 是 其 过 程 繁 
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才 程 中 像素 位 置 调整 和 图 像 选 择 。 通 过 快速 探测 宽 、 窄 视 场 图 像 的 流量 和 
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近 几 十 年 来 ， 针 对 银河 系 甚至 全 天 的 大 型 多 波段 测 光 、 光 谱 巡 天 项 目 
则 能 够 获得 发 射线 强度 ， 可 以 得 到 天 体 的 温度 、 元 素 丰 度 、 


相应 匹配 ， 


旭 像 相 减 效率 。 根 据 实验 结果 ， 减 完 连续 谱 的 图 像 中 大 概 70% 亮 星 能 被 扣除 干净 ， 并 且 能 
取 到 较 清晰 的 发 射线 天 体 结构 。 目 前 该 方法 已 用 于 搜寻 Herbig-Haro Rik (简称 HH 天 体 )， 
超新星 遗迹 候选 体 选取 等 工作 中 。 
i: 窄 波段 巡天 ; 减 连续 谱 ， 减 天 光 背 景 ， 图 像 数 据 处 理 
中 图 分 类 号 : P155 文献 标识 码 : A 


相继 开展 。 光 学 鹤 
星系 红 移 等 物理 


对 研究 激 波 发 射 源 、 超 新 星 遗迹 、H II 区 、 行 星 状 星云 等 都 有 重要 的 意义 "”。 目 前 已 有 的 罕 
波段 巡天 大 都 是 Hoc 窗 带 巡天 。 角 分 量 级 的 巡天 主要 以 SHASSA (The Southern H-Alpha 
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Sky Survey Atlas) fll VTSS" (The Virginia Tech Spectral-Line Survey) 项 目 为 代表 。 较 高 空 
间 分 辩 率 的 He 巡天 观测 有 1997 年 英 澳 天 文 台 开展 的 SHS” (Southern Halpha Survey), F 
视 场 是 5.5° x 5.5"， 其 空间 分 辨 率 可 达 1" ~ 2", 2003 年 完成 。 之 后 又 开展 了 IPHAS™ (The 


INT Photometric Ha Survey of the Northern Galactic Plane) Ha 巡天 项 


目 ， 望 远 镜 视 场 较 


小 ， 为 34 x 34， 其 空间 分 辩 率 在 1.7” 左右 ， 观 测 了 北 天 球 109 宽 的 银 盘 约 1 800° XIX. 
WHAM" (The Wisconsin Ha Mapper Survey) 完成 了 赤 续 高 于 —30° 的 3/4 天 球 的 Hoc 窄带 


覆盖 。 除 了 上 述 窜 波段 巡天 ，PTF™ (The Palomar Transient Factory, Hii 
] ) 也 完成 了 赤 纬 高 于 一 25° 的 3/4 天 球 的 Ho ite. THOSE 


洛 玛 瞬 变 源 工 


有 具体 信息 参见 表 中 


表 1 He 巡天 项 目 
项 目 名称 视 场 天 区 履 盖 范围 角 分 辨 率 Ho 中 心 波长 /A ”带宽 /A 
SHASSA 13° x 13° —90? <6<+15° 47.64" pixel ! 6 563 32 
VTSS 10° (直径 ) 6 >= —15? 1.6’pixel~! 6 563 17.5 
SHS 5.5? x 5.5? 6<2° 1" —2" pixel! 6 568 80 
IPHAS  0.3x0.3?  —5?«b«5?,290? < l< 215? 0.33"pixel ^! 6 568 95 
WHAM 1? (直径 ) 6 >= —30° l^pixel ! 6 563 20 
PTF 7.8? x 7.8? 6 >= —25? l'pixel-! 6 563 100 


Xe 5 vi FY LER RT UA RNR EEN Sz 78 RE UR RI RR RE, TERINA 
类 发 射线 天 体 候选 体 的 搜寻 ， 是 窄 波 段 图 像 数据 处 理 中 必 不 可 少 的 一 步 。 传 统 罕 波 段 数 据 
的 连续 谱 扣 除 需要 手动 去 调节 星 像 位 置 对 齐 ， 这 个 过 程 还 会 涉及 到 图 像 旋转 、 翻 转 等 问题 ， 
会 手动 选择 恒星 做 流量 校准 。 早 期 很 多 研究 甚至 是 采用 人 有 眼 进行 连续 谱 和 窒 波 段 数 据 对 比 ， 
如 Bally 等 人 在 1995 年 对 L1228 区 域 HH 双 极 喷 流 、 分 子 外 流 的 研究 ”， 以 及 Wang 等 人 
在 2003 年 对 OB1 分 子 云 中 HH 天 体 观 测 ” 都 是 采用 人 眼 对 比 ， 这 种 方式 对 于 科学 目标 的 
认证 引入 了 人 为 因素 。 之 后 减 连续 谱 工 作 ， 需 要 先 确定 好 探测 目标 区 域 ， 再 去 进行 连续 谱 扣 
除 ， 例 如 Blair 和 Long 2004 年 对 M83 超新星 遗迹 进行 光学 观测 呈 、Long 等 人 在 2019 年 对 
超新星 遗迹 的 观测 也 是 限定 在 NGC 6946 范围 内 ”。 这 种 方法 对 于 小 范围 天 区 和 目标 源 是 
可 行 的， 一 旦 涉及 巡天 大 数据 和 多 天 区 观测 目标 ， 对 目标 候选 体 的 选 定 则 比较 繁琐 ， 因 此 ， 
本 文 提出 一 种 批量 自动 化 扣除 连续 谱 的 方法 ， 相 比 传统 方法 能 够 提高 效率 。 通 过 批量 连续 谱 
的 扣除 ， 我 们 可 以 统一 从 减 完 连续 谱 图 像 中 去 搜寻 疑似 的 发 射线 天 体 。 

本 文 第 2 章 介 绍 NBS 巡天 系统 ， 包 括 观测 的 台 址 、 使 用 的 仪器 和 进展 ， 以 及 数据 的 介 


m. 第 3 章 介 绍 图 像 数据 处 理 流程 ， 主 要 介绍 连续 谱 的 扣除 ;第 4 章 进行 总 结 。 
2 NBS 巡天 简介 


NBS 是 针对 北 天 区 银 盘 开展 的 鹤 波 段 巡 天 项 目 ， 使 用 中 国 科学 院 紫 金山 天 文 台 盱 眙 观 
测 站 近 地 天 体 望 远 镜 (China Near Earth Object Survey Telescope, CNEOST)。 该 望远镜 
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采用 施 密 特 光 学 系统 ， 主 镜 口 径 120 cm， 焦 比 F/1.8: 施 密 特 改正 镜 口 径 104 cm MA 
10* x 10* CCD 相机 ， 像 元 尺寸 13.9 cm， 单 个 像素 对 应 天 空 角 1.028”"， 单 帧 覆盖 3° x 3° 视 
场 (以 下 称 为 一 个 子 天 区 )。NBS 巡天 在 每 个 农历 十 五 以 及 前 后 各 两 天 对 每 个 子 天 区 拍 5 次 ， 
每 次 曝光 180 so 

为 了 避免 大 视 场 边缘 成 像 质 量 差 ， 在 观测 中 采用 各 相 邻 子 天 区 部 分 重合 的 方法 。 如 图 吕 
所 示 ， 共 观测 9 个 子 天 区 ， 星 号 是 子 天 区 中 心 位 置 ， 灰 色 区 域 是 子 天 区 之 间 的 重 芭 部 分 ， 重 
丢 区 域 宽度 是 1?， 相 当 于 将 3° x 3? 的 视 场 均 分 了 9 块 ， 四 周 的 八 块 是 重 雪 区 域 ， 子 天 区 中 
心 外 1 一 1.5? 的 区 域 可 以 用 相 邻 子 天 区 中 心 外 0.5? —1* 的 区 域 履 盖 ， 获 得 成 像 质 量 较 高 的 
PLRK 


相对 银 纬 b/(°) 


35 25 -15 05 05 15 25 3.5 
相对 银 经 1/(°) 


RNS 9 个 子 天 区 ， 星 号 是 子 天 区 中 心 位 置 ， 灰 色 区 域 是 子 天 区 之 间 的 重 肝 部分， 重合 区 域 宽度 是 1?， 子 
天 区 中 心 外 1° —1.5? 的 区 域 用 相 邻 子 天 区 中 心 外 0.5°~1° 的 区 域 覆 盖 ， 获 得 成 像 质量 较 高 的 2° x 2” 天 区 。 


Al 子 天 区 之 间 覆 盖 示 意图 


ne 
B 


目前 NBS 系统 配备 有 [S II] 滤 光 片 、He 滤 光 


XH: m 
Ss mL am An [O m 滤 光 片 。 滤 光 片 透 过 率 曲 线 见 图 四 Hp 
[S 11] 6 T65 80 参数 见 表 2。 
He 6 563 80 由 于 本 文 使 用 的 连续 谱 数 据 为 盱眙 观测 站 波段 
a 80 图 像 ， 所 以 图 四 还 包括 盱眙 观测 站 r 波段 通 光 曲线 ， 


其 中 盱眙 观测 站 r 波段 中 心 波 长 为 6 240 Å, WA 
2500 A。NBS 观 测 的 天 区 为 北 天 区 银 盘 (—5°< b <+5°, 29°< 1 <215°) 4 1 800 平方 度 ， 从 
2016 年 开始 到 目前 为 止 ，[S IT] 波段 以 及 Ha 波段 的 第 一 遍 观 测 已 经 基本 完成 ，[O III] 波段 
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正在 拍摄 中 。 图 加 展示 了 [S Ir] 波段 数据 覆盖 图 ， 红 色 方 框 是 观测 履 盖 天 区 ， 虚 线 是 银河 系 ， 
高 银 纬 的 覆盖 是 开展 的 ENBS (extend narrow band survey). 


100 
NBS[Ba] " 
= 通 光 曲线 
E 60r NBS[S 可] 通 光 曲 线 
At 
jal 
= 40H 
NBS[O II 盱眙 观测 站 
通 光 曲线 r 波段 通 光 
20r 
0 m H 
4 000 5 000 6 000 7 000 8 000 


YA 


VE: 图 中 从 左 往 右 窄带 曲线 分 别 为 [O III], Hox, [S II] 通 光 曲线 ， 紫 色 加 粗 曲 线 为 时 眙 观测 站 r 波 段 通 光 曲 线 ， 
蓝 色 实 线 为 望远镜 反射 效率 ， 黑 色 虚 线 为 CCD 响应 曲线 。 


2 NBS 系统 滤 光 片 透 过 率 曲 线 


-120° 
| 


ik: 红色 方 框 是 观测 覆盖 天 区 ， 虚 线 是 银河 系 ， 高 银 纬 的 覆盖 是 开展 的 ENBS. 


图 3 ”NBS XIX [S IT] 波段 数据 覆盖 示意 图 


Ti 
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3 数据 处 理 流 程 


数据 处 理 流 程 包括 图 像 预 处 理 、 坐 标定 位 、 测 光 、 流 量 定 标 、 连 续 谱 扣除 、 星 表 和 罕 波 
段 数据 成 图 。 本 文 的 重点 部 分 是 流量 定 标 和 连续 谱 扣 除 。 
3.1 ”数据 预 处 理 

图 像 预 处 理 包括 减 overscan/Bias 改正 ， 除 超级 天 光平 场 ， 去 除 宇 宙 线 及 坏 像 元 等 。 预 
处 理 完 成 后 取 中 间 2° x 2° 天 区 为 目标 图 像 。 图 团 展 示 了 单 幅 图 像 的 原始 图 以 及 通过 数据 处 
理 流程 完成 预 处 理 的 图 像 展 示 。 


VE: a) [S II] 波段 原始 图 ，b) 通过 数据 处 理 流程 完成 预 处 理 的 图 像 。 


图 4 原始 图 与 预 处 理 完 的 图 像 对 比 


通过 天 文 软件 SCAMP 进行 天 体位 置 定 标 ， 使 用 UCAC-4 星 表 进行 匹配 ， 该 星 表 覆盖 
了 V AIR 波段 8~16 mag 范围 的 主要 目标 天 体 ， 天 体 总 数 多 达 113 780 093 个 ， 是 一 个 合 3 
的 全 天 高 精度 天 体 测 量 星 表 ” 。 位 置 定 标 精度 为 0.2" 左右 ， 图 卓 分 别 给 出 在 赤 经 和 赤 纬 两 
方向 的 误差 统计 ， 分 别 为 Aa = —0.02" + 0.20", Ad = —0.02" + 0.15", MARE RE ai 
3.2 ” 测 光 与 流 量 定 标 

NBS 观测 天 区 属于 密集 星 场 ， 测 光 方 法 使 用 孔径 测 光 和 点 扩散 函数 (PSF) 测 光 。 孔 径 
测 光 对 单个 目标 源 具 有 良好 的 效果 ， 但 是 当 处 于 密集 星 场 ， 就 会 对 流量 和 天 光 的 计算 产生 影 
响 。PSF 测 光 适用 于 密集 星 场 ， 这 种 情况 下 PSE 测 光 会 比 孔径 测 光 有 更 高 的 精度 。 

流量 定 标 工作 是 找到 仪器 流量 与 标准 系统 下 的 流量 的 一 个 比例 关系 ， 从 而 将 仪器 星 等 
转换 成 标准 星 等 。 由 于 在 [S IU 波段 没有 标准 星 ， 我 们 就 无 法 得 到 [S I] 波段 的 标准 流量 
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--- Aó—-0.02 +0.15 


0 05 00 05 10 1500 05 10 15 20 25 
Aa/(") 归 一 化 拟 合 


赤 纬 两 方向 天 文 位 置 定 标 的 误差 统计 ， 分 别 为 Aa = —0.02" +0.20”, Ad = —0.02" 4 0.15". 


5 ”位 置 定 标 精度 示意 图 


因此 也 不 能 求 出 比例 关系 。 这 里 我 们 提出 利用 LAMOST 光谱 数据 构建 标准 星 的 方法 : 先 改 


和 ， 从 而 得 到 
具体 步骤 如 下 : 
(1) 数据 准 


星体 被 选择 ，2 


IE LAMOST 的 仪器 流量 为 标准 流量 ， 再 用 改正 后 的 标准 流量 与 [S IT] 波段 响应 曲线 卷 积 求 


S II 波段 的 标准 流量 ， 进 而 求 得 |S II| 星 等 、 建 立 [S TI] 波段 的 标准 星星 表 。 


备 。 首 先 对 NBS 图 像 范 围 内 源 的 位 置 坐标 分 别 与 LAMOST 星 表 、IPHAS 


星 表 匹 配 ， 算 选 出 共同 目标 源 。 其 次 对 目标 源 进行 得 选 ， 1) 选择 类 型 为 恒星 的 源 ， 避 人 免 类 


) 选取 信 品 比 高 于 30 的 源 ， 保 证 源 的 精度 ，3) 选择 亮 于 15 mag 的 下 型 星 。 


最 后 保存 LAMOST 流量 数据 和 IPHAS W r 波段 星 等 数据 。 

(2) 改正 LAMOST 相对 流量 为 标准 流量 。IPHAS 巡天 观测 的 + 波段 星 等 数据 可 以 提供 
r 波段 标准 流量 ，LAMOST 仪器 流量 和 IPHAS 的 波段 响应 曲线 卷 积 求 和 获得 LAMOST 
在 r 波 段 的 相对 流量 ， 此 时 得 到 T 波 段 标准 流量 与 仪器 流量 的 定 标 常 数 ， 再 用 原来 的 仪器 流 


量 乘 上 这 个 定 标 常 数 ， 从 而 获得 LAMOST 在 r 波段 范围 内 标准 流量 。 


(3) 计算 [S I 波段 标准 流量 。 因 为 NBS 的 [S TI] 滤 光 片 波长 范围 完全 沙 在 波段 波长 


范围 内 (如 图 日 


所 示 )， 所 以 我 们 可 以 使 用 LAMOST Æ r 波段 的 标准 流量 来 计算 。[S II] 波段 
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14 000 一 一 LAMOST 相对 流量 
一 一 IPHAS r 波段 范围 
一 一 NBS [S IT] 波段 范围 
| --- IPHAS r 波段 通 光 曲线 
12 000 --- NBS [S II] 波段 通 光 曲线 
10 000* 
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E. 8 000r TRA 
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AA 
HS 蓝 色 谱 线 为 LAMOST 光谱 数据 ， 流 量 值 为 相对 流量 ， 谱 线 覆 盖 绿 色 区 域 为 IPHAS 的 波段 波长 覆盖 范 
司 ， 谱 线 履 盖 红色 区 域 为 NBS 的 [S II] 波段 波长 覆盖 范围 ， 绿 色 虚 线 为 IPHAS H r 波段 通 光 曲 线 示意 图 ， 红 
色 虚 线 为 NBS 的 [S I1] 波段 通 光 曲 线 示 意图 。 利 用 IPHAS BAH r 波段 星 等 将 LAMOST 光谱 在 v 波段 范围 


为 将 相对 流量 改 了 


E 为 标准 流量 。 


响应 曲线 通过 LAMOST 的 r 波段 标准 流量 


(4) 计算 [S TI] 


波段 星 等 


图 6 [S II] 波段 标准 星 构建 


卷 积 求 和 获得 [S IT] 波段 的 标准 流量 。 


等 ， 建 立 [S II] 波段 标准 星星 表 。 流 量 定 标 后 测 


光 精 度 如 图 加 所 


Io HFR Ta) 是 [S II] 单 幅 图 像 曝光 时 间 180 s 得 到 的 星 等 和 误差 分 布 图 ，5c 处 极限 星 等 


约 17.5 mag. Al fb) 是 


[S II] 图 像 曝光 时 间 900 s 得 到 的 星 等 和 误差 分 布 图 ，5c 处 的 极限 星 


等 可 达 18.3 mag, BRIGHT] 900 s 为 5 幅 180 s ARB. 900 s 曝光 时 间 的 星 等 基本 达到 预 
期 设计 的 5o 处 18.5 mag 的 要 求 。 


建立 [S Il] 波段 标准 星 ， 是 为 了 给 


体 的 流量 。] 


发 射线 天 体 的 流量 计算 提供 条 件 。 
像 即 为 做 完 流量 定 标 的 结果 。 对 于 发 射线 天 体 ， 
[S 2 e 
一 般 把 + WER URE E Rc BE, KETAN r 波段 数 


SS 


3.3 
减 连续 


[S IT] 波段 图 像 做 流量 定 标 ， 便 于 后 期 测量 发 射线 天 


上 述 构建 标准 星 的 方法 ， 可 以 提高 窜 波 段 流 量 定 标的 精度 ， 为 测量 减 连 续 谱 后 
对 于 做 完 流 量 定 标的 宽窄 带 图 像 ， 通 过 相 减 得 到 的 残余 儿 
直接 测 光 即 可 获得 流量 信息 。 


据 作为 连续 i eti 经 计算 盱眙 观测 站 的 r 波 段 流 量 中 [S IT] 发 射流 量 仅 占 0.25%， 所 以 从 
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ik: a) [S II] 单 幅 图 像 曝光 时 间 180 s 得 到 的 星 等 和 误差 分 布 图 ， 


曝光 时 间 900 s 得 到 的 星 等 和 误差 分 布 图 ，5c 处 的 极限 星 等 可 达 18.3 mag, 


RÆ. 


图 7 


[S I1] 罕 波 段 图 像 中 减 去 的 了 波段 


3.3.1 RADAR 8 de 


的 ， 所 以 在 做 扣除 连续 谱 时 ， 


数据 来 获取 天 光 背 景 ， 亮 源 信 


MENTATION) 图 像 ， 图 像 信 
为 0， 如 图 Bb) 所 示 。 


(2) 扣除 亮 源 信 息 。 对 输出 源 累 计 探 测 器 计数 模式 图 
值 为 1， 之 后 与 原 图 像 数据 相 乘 ， 完 成 亮 源 信 息 的 扣除 。 


图 像 中 的 [S TT] 流 


R] 


测 光 精 度 示意 图 


ta 


Beo P ER RET OG 


体 步 骤 如 下 。 


bo 处 极限 星 等 约 17.5 mag: b) [S IT] 图 像 


曝光 时 间 900 s 为 5 幅 180s 图 


可 以 忽略 不 计 。 


因为 夜 天 光 在 空间 尺度 上 存在 着 梯度 的 分 布 ， 对 于 大 视 场 巡天 望远镜 来 说 是 不 能 忽略 
背景 信息 的 扣除 。 采 用 图 像 减 亮 源 
息 的 获取 通过 Source Extract 软件 (后 文 简称 SE 软件 )， 亮 源 
是 指数 值 高 于 SE 软件 最 低 探测 阐 值 的 源 ， 最 低 探测 闵 值 由 “DETECT_THRESH” 参 数 设 置 ， 
一 般 设 置 为 3c。SE 软件 可 以 提供 源 探测 、 源 累计 探测 器 计数 、 星 像 轮 廊 的 信息 等 。 
示 了 图 像 减 天 光 背 景 过 程 ， 具 

(1) 获取 亮 源 位 置信 息 。 通 过 设置 SE 配置 文件 ， 输 昌 


3] B He 


HI ER VERUS Ht BRE X, (SEG- 


妃 只 包含 源 的 位 置信 息 即 亮 源 位 置 为 源 的 计数 ， 其 他 位 置 数据 


(3) RRAK HE Ro 


将 扣除 亮 源 信 息 后 图 


pixelx500 pixel。 天 光 背 景 如 
天 光 背 景 图 像 ， 避 免 星 云 信息 
(4) 减 天 光 背 景 。 用 原始 
Ed) 所 示 。 
NBS 图 像 覆 盖 银 盘 区 域 ， 


的 信息 的 同时 减少 噪点 。 


的 扣除 。 


图 Be) 为 下 


像 做 处 理 ， 令 


Ab EL 
Aa 


亮 源 数据 为 0， 


图 像 减 去 天 光 背 景 图 像 ， 获 得 减 完 天 光 背 景 后 的 图 像 ， 如 图 


图 像 中 大 多 包含 星云 ， 对 星云 数据 的 提取 是 一 个 难点 。 我 们 
通过 调 低 SE“DETECT THRESH” 参 数 ， 同 时 增 大 “DETECT MINAREA”， 在 得 到 更 多 
T 
FELD 


带 星云 的 [S II] 波段 图 


像 进行 中 值 平滑 处 理 。 对 于 像 NBS 这 种 大 尺 
度 图 像 ， 应 适当 增 大 平滑 尺度 以 消除 大 规模 大 幅度 的 背景 变化 ， 该 图 像 
图 Be) 所 示 ， 在 没有 星云 的 数据 中 该 方法 可 以 得 到 较为 均匀 的 


I 


F 滑 尺度 设置 为 500 


dit 


f. ABE) 为 采用 上 述 方法 后 
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g) 
a)—d) 展示 了 不 带 星云 图 
式 图 像 ，c) 天 光 背 景 图 
波段 图 像 ，f) d 


像 的 背景 提取 与 扣除 : a) AA 
像 ，d) 减 背 景 后 
zit 


的 [S IT] 波段 背景 


区 


像 ，e) 一 b) 


h) 
星云 的 S IT] RRRS b) WR ZERIT CR 
ETUR E 
图 像 ，g) 拟 合 的 背景 图 像 ，D) 减 背景 后 图 像 
图 8 减 天 光 背 景 结 果 示 意图 


取 与 扣除 : e) 含 星云 的 [S I] 
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得 到 的 背景 图 像 ， 可 以 看 到 在 右上 角 会 偏 亮 一 些 ， 这 是 因为 有 星云 信息 残留 。 如 果 用 该 图 做 
背景 扣除 ， 会 减 掉 部 分 星云 的 信息 ， 对 于 这 种 情况 ， 我 们 从 背景 图 像 中 选取 没有 星云 或 星云 
很 暗 的 区 域 (如 图 区 ) 红色 方 框 等 区 域 ) 进行 背景 拟 合 ， 图 加 中 的 g) 由 f) 红色 方 框 内 的 背景 
拟 合 而 成 。 
3.3.2 ”连续 谱 的 扣除 

为 了 提高 图 像 的 信 噪 比 ， 在 做 连续 谱 扣 除 前 ， 分 别 对 减 完 背景 后 两 个 波段 图 像 进行 中 值 
平滑 ， 这 样 可 以 使 图 像 中 激 波 天 体 结构 更 清晰 。 减 连续 谱 需 要 精确 计算 出 两 个 波段 相同 源 的 
流量 比例 因子 与 星 像 匹 配 因 子 ， 使 两 幅 图 像 源 的 流量 与 星 像 都 精确 匹配 。 流 量 比例 因子 相对 
容易 得 到 ， 星 像 匹 配 目标 是 使 两 个 波段 图 像 像 源 大 小 调 成 一 致 ， 小 像 源 通过 高 斯 滤波 函数 
G(0, a) 扩展 为 大 像 源 ， 其 中 o 高 斯 核 值 即 为 星 像 匹 配 因 子 。 因 为 星 像 的 半 高 全 宽 FWHM 
与 高 斯 函数 的 o 成 比例 关系 ， 有 具体 表示 如 下 : 


FWHM =2V/2In(2)o . (1) 


通过 调节 FWHM 可 实现 星 像 大 小 的 调整 。 连续 谱 扣除 流程 示意 图 如 图 日 所 示 ， 图 四 为 减 
连续 谱 过 程 示意 图 。 减 连续 谱 有 具体 步 又 如 下 。 


SE, python 


[s II] FERRERS 


天 光 背 景 信息 提取 


| python 


SE 


FER K/) (FWHM) 


9 ” 减 连 续 谱 流程 示意 图 


E 分 别 对 两 个 波段 数据 进行 测 光 ， 获 取 坐 标 位 置信 息 、 星 像 的 半 
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e) 


YE: a), b) 分 别 为 减 天 光 背 景 后 [S I1] 波段 图 像 和 平滑 后 的 [S II] FRA: c), d) 分 别 为 减 天 光 后 r 波段 图 


像 和 平滑 后 的 v 波段 图 像 ，e) 为 减 完 连 续 谱 的 图 像 ， 图 中 红色 框 内 均 能 够 看 到 较为 清晰 的 纤维 状 结构 ， 该 结构 在 


[S I1] 窄带 图 像 中 较为 清晰 ， 而 在 盱眙 的 v BRE 


图 10 


图 像 中 并 不 明显 。 


减 连续 谱 过 程 示意 图 
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(2) 星 像 匹 配 。 分 别 对 两 个 波段 测 得 的 PWHM 画 出 分 布 直方 图 ， 对 直方 图 进行 高 斯 拟 
合 ， 拟 合 出 峰值 (高 斯 分 布 中 风 值 ) 对 应 的 数据 ， 分 别 记 为 FWHMI 和 FWHM2， 再 通过 式 
(1) 分 别 计 算出 oy 和 os 值 ， 这 里 我 们 假设 o Kon 大 ， 由 高 斯 分 布 公式 ， 可 以 计算 星 像 匹 
ACAI o: 


e = V — (3 . (2) 


将 该 匹配 因子 作为 高 斯 卷 积 尺度 与 原 图 做 高 斯 卷 积 ， 完 成 两 幅 图 像 的 星 像 匹 配 。 

(3) 流量 匹配 。 对 匹配 好 星 像 轮 廊 的 两 幅 图 分 别 测 光 ， 相 除 得 到 各 个 源 的 流量 比例 因子 ， 
画 出 流量 比例 分 布 直方 图 ， 对 直方 图 做 高 斯 拟 合 ， 得 到 峰值 (高 斯 分 布 中 / 值 ) 对 应 的 数据 
即 为 流量 比例 因子 。 得 到 的 流量 比例 因子 通过 与 原 数据 相 乘 ， 完 成 两 幅 图 像 流量 值 的 匹配 。 

(4) 扣除 连续 谱 。 使 用 SWARP 软件 进行 两 幅 图 像 扣除 ，SWARP 软件 通过 两 幅 图 像 的 
WOS 坐标 匹配 ， 避 免 人 为 调整 像素 的 位 置 、 图 像 旋转 等 过 程 ， 处 理 速 度 快 ， 适 合 处 理 大 批 
量 数据 。 对 波段 数据 取 负 值 ， 之 后 通过 SWARP 软件 将 匹配 后 的 两 幅 图 像 求 和 ， 完 成 连续 
谱 的 扣除 。 

实验 结果 如 图 De) 所 示 ， 在 红色 框 内 均 能 够 看 到 清晰 的 纤维 状 结构 ， 该 结构 在 [S II] 
窄带 图 像 中 比较 清晰 ， 而 在 盱眙 的 r 波 段 图 像 中 并 不 明显 。 在 使 用 本 文 批 处 理 扣除 连续 谱 
后 ， 能 够 获得 较为 清晰 的 [S LI] 图 像 中 的 激 波 结构 。 我 们 分 别 对 连续 谱 扣 除 前 后 的 图 像 进 
行 测 光 ， 未 扣除 连续 谱 的 [S TT] 图 像 探测 到 63 500 个 源 ， 经 过 连续 谱 扣 除 后 的 图 像 探 测 到 
18000 多 颗 源 ， 扣 除了 原 图 像 中 70% 的 星 。 图 四 展示 了 局 部 区 域 的 星 分 布 情 况 ， 区 域 大 小 
为 2000 pixelx2000 pixel， 原 图 像 尺 寸 为 7 000 pixelx7000 pixel， 在 减 完 连续 谱 后 的 局 部 
图 像 能 够 看 到 扣除 了 大 部 分 星 。 


40 000 
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a 


iE: a) 为 连续 谱 扣除 前 局 部 区 域 的 星 像 分 布 ，b) 为 连续 谱 扣 除 后 局 部 区 域 的 星 像 分 布 。 
图 11 ”连续 谱 扣除 前 后 局 部 区 域 的 星象 分 布 


由 于 图 像 中 包含 各 种 类 型 恒星 ， 不 同类 型 的 恒星 在 匹配 宽带 和 窜 带 的 流量 时 匹配 因子 
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不 同 ， 这 就 会 导致 图 像 中 没有 发 射线 的 恒星 不 能 完全 扣除 干净 。 但 我 们 的 目标 是 从 海量 图 像 
中 ， 能 够 快速 寻找 出 含有 发 射线 天 体 的 信息 ， 一 般 具 有 延展 性 或 纤维 状 结构 。 对 于 个 别 亮 源 
言 轧 扣除 不 完全 ， 其 对 我 们 搜寻 具有 延展 结构 的 发 射线 天 体 影 响 不 是 很 大 。 目 前 该 方法 已 用 
于 搜寻 HH 天 体 和 超新星 遗迹 候选 体 选 取 等 工作 中 ， 寻 找 激 波 天 体 候选 体 的 过 程 如 图 中 所 
示 。 


o) 
He a) Jy [S TI] 波段 图 像 ，b) 为 连续 谱 图 像 ，c) 为 减 完 连 续 谱 后 [S TI] 波段 图 像 。 从 c) 可 见 疑似 激 波 天 体 。 


图 12 “寻找 激 波 天 体 候 选 体 过 程 展示 


4 总 结 


本 文 基于 窜 波 段 巡 天 项 目的 [S TI] 波段 数据 ， 提 出 一 种 通过 自动 匹配 不 同 波段 流量 数据 
与 星 像 轮廓 方法 ， 实 现 [S IT] 波段 大 视 场 图 像 批 量 扣除 连续 谱 。 我 们 引入 了 SWARP 软件 和 
Source Extract 软件 ，SWARP 软件 可 以 自动 匹配 两 幅 图 的 天 文 位 置 ， 不 需要 人 为 对 两 幅 医 
像 的 像素 与 像素 进行 对 齐 ， 避 免 了 这 个 过 程 中 像素 位 置 调整 和 图 像 旋转 ，Source Extract f 
件 能 够 快速 探测 宽 、 窜 图 像 中 的 流量 、 星 像 ， 通 过 对 其 流量 、 星 像 值 做 直方 图 统计 获得 比例 
因子 来 完成 流量 匹配 和 星 像 匹 配 。 最 后 使 用 SWARP 进行 图 像 相 减 的 批 处 理工 作 ， 相 比 传 
统 方法 能 够 提高 效率 。 根 据 实验 结果 ， 在 已 减 去 连续 谱 的 图 像 中 ， 大 概 70% 亮 星 能 被 扣除 
干净 ， 并 且 能 够 获取 到 较为 清晰 的 发 射线 天 体 结 构 。 我 们 做 连续 谱 扣 除 目 的 是 从 海量 图 像 中 
快速 寻找 出 发 射线 天 体 结构 ， 这 类 结构 一 般 有 延展 性 结构 ， 个 别 亮 源 信 息 扣 除 不 完全 对 我 们 
搜寻 具有 延展 结构 的 发 射线 天 体 影响 不 大 。 我 们 通过 第 一 步 得 选 出 含有 疑似 目标 的 候选 体 
之 后 ， 在 此 基础 上 继续 做 更 为 精细 的 图 像 处 理 ， 或 进行 发 射线 天 体 流 量 的 精确 测量 。 目 前 该 
方法 已 用 于 搜寻 HH 天 体 ， 超 新 星 遗 迹 候选 体 选取 等 工作 中 。 

目前 的 方法 还 存在 问题 ， 需 要 提高 与 完善 。 第 一 ， 对 含 大 斥 度 星云 图 像 提取 背景 时 ， 部 
分 星云 信息 残留 在 背景 中 ， 为 了 避免 将 星云 信息 扣除 ， 我 们 采用 了 小 区 域 背景 代替 完整 背景 
的 方法 ， 但 这 也 产生 背景 信息 不 准确 的 问题 ， 导 致 减 连续 谱 的 精度 不 高 ; 第 二 ， 图 像 中 包含 
各 种 类 型 恒星 ， 不 同类 型 的 恒星 在 匹配 宽带 和 窄带 的 流量 时 匹配 因子 不 同 ， 这 就 会 导致 图 像 
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中 没有 发 射线 的 恒星 不 能 完全 扣除 干净 ， 第 三 ， 对 于 一 些 亮 源 ， 星 像 匹 配 不 能 做 得 很 好 。 在 
今后 的 工作 中 ， 我 们 会 针对 以 上 问题 进一步 优化 和 完善 。 
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Abstract: Optical narrow-band observation can characterize the clear structure of the 
emission line objects, Which is important to study the star-forming region, quasar narrow- 
line region, H II region, and interstellar medium. Narrow Band Survey (NBS) is a fully 
photometric CCD survey of the Northern Galactic Plane, which targeted the inner regions 
of the Northern plane over the latitude range of —5?« b <+5° and a longitude range of 
29°< | «215? covering a total of 1800 deg?. NBS was conducted with three narrow-band 
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filters: Ha filter (80 A full width at half-maximum, central wavelength at 6563 A), [S II] filter 
(80 A full width at half-maximum, central wavelength at 6765 A) and [O III] filter (80 A 
full width at half-maximum, central wavelength at 5007 A). The continuum subtraction 
is a common method to search for emission-line objects forms the narrow-band images. 
Traditionally these have been dealt with by carefully registering the images in terms of 
position,orientation and intensity, and subtracting them from different bands images. This 
paper uses the [S II] band data of the NBS, proposes a method of subtracting continuum. 
We introduced SWARP and Source Extract, the SWARP software can automatically match 
the astronomical positions of the respective sources of the two images, avoid pixel position 
adjustment and image rotation; The Source Extract software can quickly detect the flux 
and star image in the wide-band and narrow-band image, and accomplish the flux matching 
and PSF matching by Python, finally use SWARP to do the batch processing of image 
subtraction. Compared with the traditional method, this method can improve the efficiency. 
At present, this method has been applied in searching for the Herbig-Haro objects and 


Supernova Remnant candidates. 


Key words: survey in narrow band; continuum subtraction; sky-background subtraction; 


image processing 


